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We give a systematic treatment of a spin 1/2 particle in a combined electromagnetic field and a 
weak gravitational field that is produced by a slowly moving matter source. This paper continues 
previous work on a spin zero particle, but it is largely self-contained and may serve as an introduction 
to spinors in a Riemann space. The analysis is based on the Dirac equation expressed in generally 
covariant form and coupled minimally to the electromagnetic field. The restriction to a slowly moving 
matter source, such as the earth, allows us to describe the gravitational field by a gravitoelectric 
(Newtonian) potential and a gravitomagnetic (frame-dragging) vector potential, the existence of 
which has recently been experimentally verified. Our main interest is the coupling of the orbital 
and spin angular momenta of the particle to the gravitomagnetic field. Specifically we calculate 
the gravitational gyromagnetic ratio as g = 1 ; this is to be compared with the electromagnetic 
gyromagnetic ratio of g e = 2 for a Dirac electron. 
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1. INTRODUCTION 

Classical systems in external gravitational fields have 
been studied for centuries, and recently the existence of 
the gravitomagnetic (or frame-dragging) field caused by 
the earth's rotation has been observed by the Gravity 
Probe B (GPB) satellite [TH1]. GPB verified the predic- 
tion of general relativity for the gravitomagnetic preces- 
sion of a gyroscope in earth orbit (42 mas/yr) to better 
than 20% [5]. Previously, observations of the LAGEOS 
satellites also indicated the existence of the gravitomag- 
netic interaction via its effect on the satellite orbits |5, 6J. 
Analysis of the LAGEOS data involves modeling classi- 
cal effects to very high accuracy in order to extract the 
gravitomagnetic effect, and the accuracy of the results 
has been questioned by some authors [7J. Analysis of 
the GPB data also requires highly accurate modeling of 
classical effects [5]. 

While gravitomagnetic effects are generally quite small 
in the solar system it is widely believed that they may 
play a large role in jets from active galactic nuclei, so 
their experimental verification is of more than theoretical 
interest [5]. 

At the other end of the interest spectrum extensive the- 
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oretical work has been done on quantum fields in classical 
background spaces, the most well known being related to 
Hawking radiation from black holes [TJ r5HTTj . However 
it is important to keep in mind that Hawking radiation 
has never been observed. 

Interesting experimental work has also been done on 
quantum systems in the earth's gravitational field, such 
as neutrons interacting with the earths Newtonian field 
and atom interferometer experiments aimed at accurately 
testing the equivalence principle and other subtle gen- 
eral relativistic effects [T2HT4] . There has been some dis- 
cussion of attempts to see gravitomagnetic effects with 
these devices but such experiments would be quite diffi- 
cult due to the small size of the effects and the similarity 
to classical effects of rotation; this is to be expected since 
gravitomagnetism manifests itself in a way that is quite 
similar to rotation, hence the appellation "frame drag- 
ging." Laboratory detection of gravitomagnetic effects 
on a quantum system would clearly be of fundamental 
interest. 

In this work we give a systematic treatment of a spin 
1 /2 particle in a combined electromagnetic field and weak 
gravitational field; this continues the work of reference 
|15j . We describe the particle with the generally covari- 
ant Dirac equation in a Riemann space, minimally cou- 
pled to the electromagnetic field in the standard gauge 
invariant way |16l 117) . The weak gravitational field is 
naturally treated according to linearized general relativ- 
ity theory, and we also assume a slowly moving matter 
source, such as the earth [TS1120] . Within this approxi- 
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mation the gravitational field is described by a gravito- 
electric (or Newtonian) potential and a gravitomagnetic 
(or frame-dragging) vector potential, and the field equa- 
tions are quite analogous to those of classical electromag- 
netism. We thus refer to it as the gravitoelectromagnetic 
(GEM) approximation. Our special emphasis throughout 
this paper is on the gravitomagnetic interaction. 

The paper is organized as follows. After brief review 
comments on the GEM approximation (section 2) and the 
Dirac equation in flat space (section 3) we give a detailed 
discussion of generally covariant spinor theory and the 
Dirac equation, using the standard approach based on 
tetrads (sections 4 and 5). We then obtain the limit of 
the Dirac Lagrangian and the Dirac equation for a weak 
gravitational field and discuss its interpretation in terms 
of an energy- momentum tensor (section 6). 

Our discussion of generally covariant spinors and 
the generally covariant Dirac equation is largely self- 
contained, and may serve as an introduction to the sub- 
ject for uninitiated readers. In section 6 we also observe 
that the non-geometric or "flat space gravity" approach 
of Feynman, Weinberg and others does not appear to be 
completely equivalent to linearized general relativity the- 
ory in its coupling to spin |5T]. We have not found this 
discussed in the literature. 

Using the weak gravitational field results we then ob- 
tain the non-relativistic limit of the theory (section 7). 
We do this by integrating the interaction Lagrangian to 
obtain the interaction energy of the spinor particle with 
the electromagnetic and the GEM fields, and from that 
obtain the non-relativistic interaction energies. This al- 
lows us to read off, in a simple and intuitive way, the 
interaction terms that one could use in a non-relativistic 
Hamiltonian treatment. In particular we obtain (section 
8) the usual anomalous g-factor of the electron g e = 2 and 
the analogous result for the gravitomagnetic g-factor of 
a spinor, which is g = 1. 

Section 8 also contains brief comments on the numeri- 
cal value of some interesting and conceivably observable 
quantities such as the precession of a spinning particle 
in the earth's gravitomagnetic field and its relation to 
the precession of a macroscopic gyroscope; such preces- 
sion appears to be universal for bodies with angular mo- 
mentum. The phase shift in an atom interferometer is 
also mentioned as an experiment that could, in principle, 
show the existence of the gravitomagnetic field. 

Lastly it is worth noting what we do not do in this 
paper. We study the effect of the gravitational field on 
a quantum mechanical spinor but not the effect of the 
spinor on the gravitational field; thus the work does not 
relate to quantum gravity or quantum spacetime |22j . 
Similarly we do not consider torsion, in which the affine 
connections have an anti-symmetric part and are not 
equal to the Christoffel symbols. Torsion does not prove 
necessary in our discussion, but some authors believe it 
is necessary in describing the effects of spin on gravity 



2. THE GRAVITOELECTROMAGNETIC (GEM) 
APPROXIMATION 

In previous work we discussed linearized general rela- 
tivity theory for slowly moving matter sources like the 
earth [T3 EH [20] • Here we summarize the results very 
briefly. The metric may be written as the Lorentz metric 
plus a small perturbation, 

We use coordinate freedom to impose the Lorentz gauge 
condition 

(V - \^M V = 0, (2.2) 

where the single slash denotes an ordinary derivative. 
Then the field equations of general relativity tell us that 
the metric perturbation may be written as 

/ 20 h 1 h 2 h 3 \ 
, _ | h 1 20 1 , _ _ , , _ , fc 
M!/ — I h 2 20 ' — — ' 

\ h 3 2(f) J 

(2.3) 

where is the Newtonian or gravitoelectric potential and 
h k f> h is the gravitomagnetic potential. For slowly mov- 
ing sources the field equations and the Lorentz condition 
become 

V 2 = 4ttGp, V 2 hP = -16irGpv j , 40 - V • h = 0, h = 0, 

(2.4) 

where p is the source mass-energy density and v J is its 
velocity. 

The physical fields, which exert forces on particles, are 
the gravitoelectric (or Newtonian) field and the gravit- 
omagnetic (or frame-dragging) field, which are defined 
by 

g = V0, n = Vxft. (2.5) 

We call this equation system the gravitoelectromagnetic 
or GEM limit because of its close similarity to classical 
electromagnetism. 



3. FLAT SPACE DIRAC EQUATION AND THE 
NON-RELATIVISTIC LIMIT 

In this section we discuss the Dirac equation in the 
flat space of special relativity and recast it into a 
Schroedinger equation form (SEF), which provides one 
convenient way to obtain the non-relativistic limit |16j . 
The SEF is exact and involves only the upper two com- 
ponents of the spinor wave function - the relevant compo- 
nents for positive energy solutions in the non-relativistic 
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limit. One reason for doing this is to serve as a ba- 
sis of comparison for the alternative method we will use 
in section 6 when we discuss gravitational interactions. 
Throughout this section 7^ denotes the flat space Dirac 
matrices [TT)ll2"4"] . 

The Dirac Lagrangian and the Euler-Lagrange equa- 
tions that follow from it are 

L = atpii^dfi — m)4> + btp(—i^ d ll — m)ip — eA^-f^ip, 

(3.1a) 

(3.1b) 

The spinor and its adjoint are considered independent in 
obtaining (3.1b). The constants a and b are arbitrary, 
so long as a + b 7^ . The 7^ obey the flat space Dirac 
algebra, 



{7", 7°} = 2,^1. 



(3.1c) 



The adjoint spinor is assumed to be related to the spinor 
by a linear metric relation, ip — ip'M where M is to be 
determined; consistency of the equations (3.1b) is then 
assured if M obeys 

M~ 1 ^ M = Y, M- 1 =M = j°, $ = ^7°- (3.2) 

Eq. (3.2) is easy to verify for the choice of gamma ma- 
trices given below in (3.4). 

The Hamiltonian form of the Dirac equation will be 
useful for studying interaction energies in this section. It 
is gotten by multiplying (3.1) by 7 to obtain 

id t ip = fimip + V + a -Hip, /3 = 7°, a = 7°7 fe , p = -iV 

(3.3) 

Pauli's choice of gamma matrices is natural for our later 
discussion of the non-relativistic limit, 



/3 = 7 ° = 
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(3.4) 

Next we break the 4-component wave function ip into 
two 2-component Pauli spinor wave functions and also 
factor out the time dependence due to the rest mass by 
substituting 



which leads to the coupled equations, 



(3.5) 



id t * = + {a ■ 11)93, id t (j> + 2mip -Vip= (a ■ Pi)*. 

(3.6) 

We are interested in * so we solve for (p , and obtain 
symbolically, 

id t V = + (a ■ n)(2m - V + id t )- l (a ■ n)*, (3.7a) 
ip = (2m- V + idt)- 1 ^- fl)*. (3.7b) 



The inverse operator (2m — V + idt)~ 1 may be defined by 
its expansion in the time derivative, as discussed in Ap- 
pendix A. Note that Eq. (3.7a) is an exact equation for 
*, although it is of infinite order in the time derivative. 

For the special case of a free particle the operator fac- 
tors on the right side of (3.7a) commute and it becomes 
simply 



idtip = (i&t + 2m) 1 fPip. 



(3.8) 



However the operators on the right side of (3.7a) will not 
in general commute unless the field is constant. 

In a low velocity system the time variations of * and 
V are associated with non-relativistic energies, which are 
much less than the rest energy m, so we may approximate 
(3.7a) by 



idtV = Vip + ^ - S qr, 
2m 



(3.9) 



This is the Schroedinger equation for spin 1/2 parti- 
cles, often called the Pauli equation. The Pauli equation 
shows clearly how the spin and orbital angular momen- 
tum interact with the magnetic field. Pauli spin matrix 
algebra leads to an illuminating form for (3.9): to lowest 
order in e, 



n 2 

2m 
f_ 
2m 



eB ■ a 

2m 
eA ■ p 



9- — *, 

Im 



(3.10) 



where we have used the Lorentz gauge in which V • A = 
—A and assumed the Coulombic A° has negligible time 
dependence. The gyromagnetic ratio or g-factor of a par- 
ticle or system is defined in terms of its magnetic moment 
ft and angular momentum J by ft — g e (e/2m) J ; thus, 
from (3.10), the fact that the energy is —p, ■ B, and the 
electron spin of S = a/2 it is evident that the electron 
g-factor is g e = 2. 

The relative coupling of the spin and orbital magnetic 
moments is made most clear if we consider a magnetic 
field that is approximately constant over the size of the 
system, in which case we can choose A = (B x r)/2 and 
find from (3.10) 



-n 
V 



eB 



2m 2m 
£=0-/2, L = fxp. 



(2S + L)% 



(3.11) 
1 for the 



That is g e = 2 for the electron spin and g, 
orbital angular momentum. 

Equation (3.7a) may be expanded to higher order to 
study such things as hyperfine structure and relativistic 
corrections in the hydrogen atom spectrum |25j . That is 



(ff-n)s 

2m 



{a-Il)(id t -V)(a-Il) 



4m 2 



(3.12) 
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However an important problem and caveat is that the 
wave function \fr in (3.12) is only the upper half of the 
Dirac wave function, so the quantity that must be nor- 
malized is |^| 2 + |ty?| 2 rather than \ip s \ 2 for a Schroedinger 
or Pauli wave function ip s . Thus to insure Hermiticity 
and conserve probability one must renormalize the wave 
function as discussed in detail in ref. [25j . It is for this 
reason that we will adopt an alternative and conceptually 
simpler approach to the non-relativistic limit in section 
7. 



vector we define a rule for transplanting a spinor from x 
to a nearby point x + dx , 



ip*{x + dx) = ip(x) — T fl ip(x)dx ,J 



(4.5) 



The matrices T a are variously called spin connections, 
affine spin connections, or Fock-Ivanenko coefficients. 
The covariant derivative is then defined in terms of the 
difference between the value of the spinor and the value 
it would have if transplanted to the nearby point. That 
is 



4. GENERALLY COVARIANT SPINOR 
THEORY 

The gravitational interaction of a spinor may be ob- 
tained most easily by making the Dirac Lagrangian (3.1a) 
and Dirac equation (3.1b) generally covariant. To do this 
we adopt the standard approach of using a tetrad of basis 
vectors in order to relate the generally covariant theory 
to the special relativistic theory in Lorentz coordinates 
[T4"l I17j. This is a most natural, almost inevitable, ap- 
proach since Dirac spinors transform by the lowest di- 
mensional representation S of the Lorentz group; that is 
tp' = Sip . 

Two properties of the Dirac Lagrangian and Dirac 
equation must be modified to obtain a generally covari- 
ant theory: the Dirac algebra in (3.2) must be made 
covariant and the derivative of the spinor in (3.1) must 
be made into a covariant derivative. We will discuss both 
in detail. 

The Dirac algebra (3.1c) is easily made covariant by re- 
placing the Lorentz metric rf v by the Riemannian metric 



{ 7 M ,7 Q } = 2g^7. 



(4.1) 



A set of 7 P matrices that satisfy (4.1) is easily con- 
structed by using a set of constant j b that satisfies the 
special relativistic relation (3.2) and a tetrad field nor- 
malized by the usual tetrad relations 



^b c &bfj.v — 'hb, g — e c e d'/ ■ 



(4.2) 



Here the Greek indices label components of the tetrad 
vectors and Latin indices label the vectors. In terms of 
a convenient set of constant Dirac matrices 7 , such as 
those in (3.4), we define the 7^ by 



(4.3) 



It then follows from (3.1c) and (4.2) that the 7 M satisfy 

Y} = «{ 7 b , n = <W>7 = 2gTL (4.4) 

The covariant derivative of a spinor is defined so as to 
transform as a vector under general coordinate transfor- 
mations and as a spinor under Lorentz transformation of 
the tetrad basis. As with the covariant derivative of a 



ip{x)\\,jdx u = [4>(x) + ip(x)\ v dx v ] — [ip(x) — T u {x)ip{x)dx v } 
= [ip(x)\ v + T v {x)ip{x)]dx v , 
i>\\ v = ^\ v + T v ^={d v +T v )^ = D u i). (4.6) 

Here the double slash denotes a covariant derivative. 
Since the spinor covariant derivative must transform as a 
vector under coordinate transformations and as a spinor 
under Lorentz transformations of the tetrad basis, we 
have 



dx v 
dx'f* 



(4.7) 



It follows from (4.6) and (4.7) that the spin connections 
must transform according to 



r' 



dx v 
dx'f* 



(4.8) 



The transformation (4.8) is formally similar to that of the 
affine connections used for vector covariant derivatives. 

The covariant derivative of an adjoint spinor follows 
easily from that of a spinor in (4.6); we ask that the 
inner product ip\ of a spinor \ an d an adjoint spinor tp 
be a scalar and thus have a covariant derivative (V'x)llAi 
equal to the ordinary derivative (tpx)\^j an( i we also ask 
that the product rule hold for both the ordinary and the 
covariant derivatives. The result is 



(4.9) 



The same idea leads to the covariant derivative of a 
gamma matrix, with only a bit more algebra; that is we 
ask that the expression ('(/'7 M x)||a be a second rank ten- 
sor and that it obey the product rule of differentiation, 
and find from (4.6) and (4.9) 



^ - V 4- 1 fl 



(4.10) 



This expression plays an important role in obtaining the 
spin connections in the next section. 



5. COVARIANT DIRAC LAGRANGIAN AND 
DIRAC EQUATION 

In this section we give a covariant Lagrangian and ob- 
tain the covariant Dirac equation. In the process we get 
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a relation between the spinor and its adjoint (i.e. a spin 
metric) and evaluate the spin connections. 

The choice of a covariant Dirac Lagrangian L, and its 
associated Lagrangian density £, is rather obvious from 
the flat space Lagrangian in (3.1), 

L = a^^VllA* - «#) + K-m^ImV 1 - i>m)ip, C = i/gL. 

(5.1) 

Coupling to the electromagnetic field will be included 
later. The 7 M denotes the covariant Dirac matrices (4.3) 
throughout this section. The Dirac equations for the 
spinor and the adjoint spinor follow directly as the Euler- 
Lagrange equations of the Lagrangian density C with ip 
and if> treated as independent variables, 

(a + - mip) + iby^ = (5.2a) 

(a + 6X^11^ + mij)) + ia^rf^ = 0. (5.2b) 

For simplicity we assume that the spin connections, un- 
specified up to this point, may be chosen so that the di- 
vergence of 7° that appears in (5.2) vanishes, 7^, = 0. 
The covariant Dirac equation is then the obvious gen- 
eralization of the flat space equations (3.1). The spin 
connections will be obtained below. Also for simplicity 
and symmetry we choose henceforth a = b = 1/2; this 
will prove convenient later. 

Next, as in flat space in section 3, we ask that there 
be a relation between the adjoint and the spinor, -tp = 
ip^M, such that the two equations (5.2) are consistent. 
Manipulating (5.2a) we get for the adjoint, 

~ i'ipM^M^ - itpf^ - ipm = 0, 

7^ sM^V 1 M, = M~ 1 T / lM. (5.3) 

We then compare (5.3) with (5.2b), written as 

-i^y + *^r„y - $m = o, (5.4) 

and see that M must satisfy the following two equations 

7* i = y = M-V t M, (5.5a) 
-r„ = f „ = M^T^M + M-^M. (5.5b) 

Eq. (5.5a) may be written in terms of flat space j b as 

e^ b = e£M- 1 7 6t M. (5.6) 

Thus it is obvious that we should ask j b = M'^M, 
which is the same as in the case of flat space and special 
relativity (3.2), so we may choose M _1 = M = 7 . Then 
the derivative of M is zero, and it is easy to verify that 
the choice M -1 = M — 7 also satisfies (5.5b). 

Our remaining task is to obtain specific spin connec- 
tions r a . To do this we make the natural demand that 
7^ have a null covariant derivative, so from (4.10) 

7ff^7r Q + {^}/ + [r Q ,7-]=0. (5.7) 



This guarantees that the divergence vanishes, 7^, = 0, 
as we have already mentioned. However it is a stronger 
demand analogous to the standard demand in general rel- 
ativity that the metric have a null covariant derivative, 
which forces the affinc connections to be the Christoffcl 
symbols. Note also that T a is obviously arbitrary up to 
a multiple of the identity, which we will suppress hence- 
forth. 

To solve (5.7) we express 7^ in terms of flat space gam- 
mas 7 b as in (4.3) and rewrite (5.7) as 

<|| Q 7 6 + [r Q ,7"]< = 0. (5.8) 
Multiplying this by the inverse tetrad matrix we get 

[F Q , 7 c ] = -^< ||ct 7 b - (5-9) 

We next note the well-known commutation relation 
on the sigma matrices, which are defined as <7 ah = 

[<7 a6 , 7 C ] = 2i{^jf c - 7 V C ), (5.10) 

From (5.10) it is evident that we should seek a solution 
that is proportional to a times a product of the tetrad 
and its derivatives. It is easy to verify that the specific 
choice 

= \^<\\J & \ (5-ii) 

satisfies (5.9) and thus serves as the spin connection. 

We thus have obtained a generally covariant theory in 
which the Lagrangian, the Dirac equations, the relation 
of the spinor to its adjoint, and the spin connections are 
generally covariant and consistent. 

Finally we include coupling to the electromagnetic field 
in the usual minimal coupling way, that is by substitut- 
ing iD^ — > iD^ — eA^; this gives the complete covariant 
Lagrangian 

L = ^(h^V'Hm - mip) + ^(-#11^ - i' m )' l P - eA^^ip. 

(5.12) 

We will study the weak gravitational field limit of this in 
the next section. 



6. LINEARIZED THEORY FOR WEAK 
GRAVITY 

In this section we use the results of section 5 for covari- 
ant spinor theory to work out the weak field linearized 
theory. This is done by setting up an appropriate tetrad 
and using it to expand the Lagrangian (5.12) to lowest 
order in the metric perturbation. The result is that there 
are three interaction terms in the Lagrangian, one asso- 
ciated with the spin coefficients and the second with the 
alteration in the 7 M caused by gravity. Remarkably the 
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first vanishes in the linearized theory, while the second 
corresponds to an interaction via the energy momentum 
tensor, as intuition should suggest. The third term is a 
cross term between the weak gravity and electromagnetic 
fields. 

In a space with a nearly Lorentz metric (2.1) it is natu- 
ral to choose a tetrad that lies nearly along the coordinate 
axes, 

= # + el = 6l-wi, (6.1) 

where w% is a small quantity to be determined. From 
the fundamental tetrad relation (4.2) it follows that we 
should choose a symmetric w^ v — — (l/2)h /M/ and thus 
have a tetrad and 7^ matrices given by 

r = [k - (1/2) W = T ~ (l/2)<7°- (6.2) 

Since Greek tensor indices and Latin tetrad indices are 
intimately mixed in the linearized theory we will not dis- 
tinguish between them in this section. 

To evaluate the spin connections (5.11) with the tetrad 
(6.2) we need the Christoffel symbols and the covariant 
derivative of the tetrad to first order in h^ u , 

{ A l} = ( 1 /2)(^ Ul + V,.-^) I 

<||^ = (l/2)(V|a-V |l/ )- ( 6 - 3 ) 

From (5.11), (6.2) and (6.3) we obtain the spin connec- 
tions, 

= \^<\\^ b = -h^. (6.4) 
Thus the Dirac Lagrangian (5.12) becomes, 

L = ^(i7 M Vv ~ + ^("^ImT^ 1 - i>rn)i> - eA^^ip 

(6.5) 

with given in (6.4). The first line is the Dirac La- 
grangian in fiat space (3.1a), and the other three terms 
are gravitational interactions that we now address. 

The first interaction term in the second line of 
(6.5), due to the spin connections, contains the anti- 
commutator {7 M ,r M }. With the use of the symmetry 
of /i M „, the Dirac algebra (3.1c), and the operator iden- 
tity [AB, C] = A{B, C} - {A, C}B it is straightforward 
to verify the following two expressions, 

K b \ a a ab r = i(h {b - h a b]a )f, (6.6) 
and thereby see that 

{r,r,}= l l h, bla {r i a ab } = 0- (6.7) 



Thus the interaction term containing the spin connec- 
tions in (6.5) vanishes, which is a remarkable simplifi- 
cation. It should be stressed that this is only true to 
first order, and the spin connections will generally be of 
interest in the full theory. 

There remains in the Lagrangian (6.5) only interac- 
tions due to the modification of the by gravity in 
(6.2); L may now be written as 

L = ~V>(£7Vv - mil)) + -(-ii>\fj,1 - V>m)V> - eA^^ip 

(6.8) 

The quantity in brackets in (6.8) is the appropriately 
symmetrized energy-momentum tensor T 1 ^ for the Dirac 
field interacting with the electromagnetic field; that 
is, the gravitational interaction Lagrangian may be ex- 
pressed as 

Lig = -\h» a [\^TW\, ~ eA^) - i(«/v + eA^TiA 

(6.9) 

The energy momentum tensor is discussed further in Ap- 
pendix B. 

The interaction (6.9) consists of the inner product of 
the field h^ v with the conserved energy-momentum ten- 
sor T^ v ; this coupling is in close analogy with the elec- 
tromagnetic coupling between the field A^ and the con- 
served current = etp^ip m (6.8). Feynman has em- 
phasized this analogy and developed a complete "flat 
space" gravitational theory, with gravity treated as an 
"ordinary" two index (spin 2) field and formulated by 
analogy with electromagnetism, at least to lowest order 
|21j . The geometric interpretation of gravity is thereby 
suppressed or ignored. Weinberg has similarly stressed 
that the geometric interpretation of gravity is not es- 
sential [T5J |5T] . Schwinger also has used a similar and 
probably equivalent non-geometric methodology called 
source theory to obtain the standard results of general 
relativity theory, including the precession of a gyroscope 
due to the gravitomagnetic field [26] . However there is a 
problem with relating the geometric and non-geometric 
viewpoints, in that the Euler-Lagrange field equations 
are based on the Lagrangian density C^/gL = (1 + h/2)L 
and not the Lagrangian L, so there is an additional in- 
teraction term (h/2)L in the geometric theory that is not 
present in the non-geometric theory; the equivalence of 
the Feynman approach to the linearized geometric ap- 
proach is thus spoiled whenever the additional term does 
not vanish. 

The difference between the Dirac equation per our ge- 
ometric development and that which one would obtain 
from the non-geometric approach is easy to see. The 
Dirac equation that follows from (6.8) with C = yfgL = 
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(l + h/2)L is 
7 M (#|/x - eA^) - rmp 

(6.10) 

The last term on the right containing h\ v would not be 
present in the non-geometric approach. This will be dis- 
cussed lurther in section 7. 

In summary of this section, the Lagrangian (6.8) con- 
tains the interaction of the Dirac field with the electro- 
magnetic field to all orders and the interaction with the 
gravitational field only to lowest order; (6.10) is the cor- 
responding Dirac equation. We will discuss the interac- 
tion energies further in the following section in which we 
consider the non-relativistic or low velocity limit of the 
theory. 



7. NON-RELATIVISTIC LIMIT 

We wish to use the results of the previous sections to 
obtain a non-relativistic limit of the theory and calcu- 
late in a simple way some interesting properties of a spin 
1/2 particle such as the electromagnetic g- factor and its 
gravitomagnetic analogue. The most familiar approach 
to this problem is to work with the upper two compo- 
nents of the Dirac wave function as we did in section 3, 
and take the non-relativistic limit [THl US] . However the 
alternative approach we use in this section is conceptu- 
ally simpler and avoids the problems of rcnormalization 
and Hermiticity that occur in the approach of sec. 3. The 
basic idea is to integrate the interaction Lagrangian over 
3-space to get the interaction energy, then put the en- 
ergy expression with Dirac 4-spinor wave functions, into 
a form using Pauli 2-spinor wave functions, all in the low 
velocity limit. 

In this section we will always work in nearly flat space 
with Lorentz coordinates; the Dirac 7 P will be those of 
flat space and no hat will be used. Moreover for simplicity 
we will work in the Lorentz gauge for both the electro- 
magnetic and GEM fields, and take both the Coulomb 
potential ^4° and the Newtonian potential 4> to have neg- 
ligible time dependence; that is A — — V • A — and 
40 = V • h = 0. This is quite appropriate, for exam- 
ple, for electromagnetic interactions in atoms and GEM 
interactions on the earth. 

To illustrate the method we first consider only the elec- 
tromagnetic interaction in flat space; the results will be 
the same as those in section 3, in particular g e = 2. The 
interaction Lagrangian and the interaction energy are, 
from (6.8), 

Ljem = -eA^VVV) = -A,j*, (7.1a) 
AE EM = - [ L IEM d 3 x. (7.1b) 



For the Dirac ip we use a convenient device, an expansion 
in terms of free positive energy Dirac wave functions on 
the mass shell. That is 

^ = E / S^fkP,s)V P "x a u{p,s)], 

s=l,2"' ^ ' 

E 2 = (p°) 2 =p 2 + m 2 . (7.2) 

The positive energy wave functions do not form a com- 
plete set, but the approximation (7.2) should be quite 
good for distances much larger than the Compton wave- 
length, h/m; (7.2) is our fundamental assumption. A key 
idea in the calculation is to express the Dirac 4-spinor 
u(p, s) in terms of a Pauli 2-spinor \s [2"T] . 

e-^uip, s) = e~^ a J ( l p )xs- (7.3) 

Correspondingly we express the non-relativistic Pauli 
wave function as 

*= 5 a /^ /(p '^ e *""" x - (7 - 4) 

In terms of the above expressions (7.2) and (7.3) the in- 
teraction energy (7.1b) is 

[e J d 3 xe tip '-- p " )xa u(p' : s') 1 ' i u(p,s)A fl }. 

(7.5) 

The bracket in (7.5) corresponds to scattering of a free 
Dirac spinor by an external field, which is equivalent 
to scattering by an infinitely heavy source particle. It 
contains all the information about the spin interaction 
and corresponds to the diagram in fig. 7.1: the parti- 
cle leaves the wave function blob with 3-momentum p, 
scatters from the external field via the QED vertex am- 
plitude into momentum p' , and then reenters the wave 
function blob. The electron remains on the mass shell, 
corresponding to zero energy transfer, which is consistent 
with a non-relativistic wave function. We denote the 4- 
momentum transfer by g M = p'^ — p^, with qo = . The 
magnitude of the allowed 3-momcntum transfer q is lim- 
ited by the width of the function f(p, s) in momentum 
space. 

It is now straightforward to calculate the bracket in 
(7.5). We split it into 2 parts, /i = for the electric 
interaction and /it = j for the magnetic interaction. For 



have 




f(p',s') 



d 3 xAjt 



u{p , s )j J u(p, s) 



I, 



E + m 

2m 
a ■ pi 
E + m 



X' s > 

er J 
a J 



/ 

B-p 
E+m 



j d 3 xe^ x ° (£-^ x \,[2p- A + q- A + iqX A-a] Xs 

(7.s; 



FIG. 1: The electron in the wave function scatters from the 
field and back into the wave function. 



the electric part we have 



e / d 3 xA e iq - x u{p',s')^u{p,s) 



d'xe^Ao'^^^' 



2m 



a ■ $ 



E + M 



I 

B-p 



= e / d 3 xe iq * x 



E q ■ p iqxp-a 

to 2m(E + m) 2m(E + m) 



Xs- 

(7.6) 



The first term in the bracket in (7.6) is the obvious charge 
coupling to the Coulomb field. The second and third 
terms may be simplified. First, because there is no energy 

transferred p = p' , from which it follows that p ■ q = 
— q 2, /2. Secondly the vector q multiplying the exponential 
may be replaced by iV operating on the exponential, 
after which integration by parts allows us to replace it 
with — iV operating on the function A$; that is we may 
replace qAo — > — iV Aq. Thus the second term vanishes 
since V 2 ^4o = in a charge free region for the Lorentz 
gauge. What remains is, to order 1/m 2 , 



< l d 3 xA e lq <* x °'u{p\s') 1 u{p,s) 

d 3 e^ x " [e{ X \,Xs) + ^V0 C x p- {xUxs) 

(7.7) 

The second term in (7.7) is clearly a nonlocal fine struc- 
ture correction, which we mentioned in sec. 4 and which 
will not concern us further in the present work [25] . 

The fi = j magnetic part of the interaction (7.5) is 
handled in exactly the same way as the electric part. We 



We then replace q — > — iV as discussed above and see 
that the second term in the bracket vanishes in a gauge 
with V • A = 0, and we are left with 



d 3 xAje iqc " x " u(p' ', s')-f j u(p, s) 

-J d 3 xe^ xa [£-^xlW-A + VxA-a] Xs 



d 3 xe tq " x ° 



^P-A(xU) + -^B-(xUxs) 



(7.9) 

Finally we combine (7.7) and (7.9) and substitute into 
(7.5) to obtain, to order 1/m, 



s,s' = l,2 



d 3 p d 3 p' 
(2tt) 3 (2tt) 3 



^.ce-^'-^xlleAo - -p- A—B ■ a] Xs 

= { d 3 x¥\eA ~ —p- A~ -^—B • (7.10) 
J m 2m 

This is the same result that we discussed in section 3, so 
we have thus verified that our present approach repro- 
duces the usual result for the electron g factor, g e = 2. 

We now work out the non-relativistic limit of the grav- 
itational interaction in (6.8), following the same proce- 
dure as for the electromagnetic interaction; we will not 
include the product of the electromagnetic and gravita- 
tional fields, that is the cross term in (6.8). The algebra 
is a bit lengthier but equally straightforward. As with 
the Lagrangian and energy for the electromagnetic case 
in (7.1) we have for the gravitational case 



f*(p',s')f(p,s) 



Lig = -\K^\ AE g = - / L IG d 3 x 



(7.11) 



where T^ v is given in (6.9). It is convenient to write T^ v 
in close analogy with the electromagnetic current, as 



(7.12) 
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Note the relation between the electromagnetic and the 
gravitational interactions, 



A^h^/2, yo7"(£tf"). 



(7.13) 



Then Ai?c is, in analogy with (7.5), 

d 3 p d 3 p' 
(2tt) 3 (2tt) 

s *&4-P-)x' , fi ( J / ) S ') 7 M(^ + g72)u(p, S )(V/2)] 

(7.14) 



d 3 G 



As with the electromagnetism calculation we split the 
gravitational interaction into two parts, the gravitoelec- 
tric for hoo = ha = 2<p and the gravitomagnetic for 
liQj = hjQ = h? . The gravitoelectric part of the bracket 
in (7.14) involves the same spin products as encountered 
with the electromagnetic calculation in (7.7) and (7.9), 
and after some algebra we obtain, to order l/m 2 , 



J d 3 xe i{p '- 



d 3 xe lip '« 



p «> x u{p', s'){^E + {p> + \W}u{ Pl s) 





f E 2 

— <5 
y m 


> + A E 2* 

4m z 


+( 


m 


2m 


mxl 




m z 



4m 2 



(7.15) 



A word is in order about the physical interpretation of 
the gravitoelectric result (7.15). The term m<f> is of course 
the expected Newtonian energy; the factor (1 + 2p 2 /m 2 ) 
occurs also in the analysis of a spin zero system in ref. 
[15j . and is approximately the Lorcntz transformation 
factor between the potential in the lab frame and the 
moving frame of the particle; thus (1 + 2p 2 /m 2 )</> is the 
Newtonian potential seen by the moving particle. The 
last term in the bracket has the same form and is the 
gravitational analog of the fine structure term in the elec- 
tromagnetic energy (7.7), except of course for the differ- 
ent coefficient. We will not be concerned further with 
the higher order terms in (7.15) and will henceforth keep 
only the lowest order term <f> in the bracket. 

We turn finally to the gravitomagnetic part of the in- 
teraction (7.14), which is our main interest in this work. 
The gravitomagnetic part of the bracket, proportional to 
h> , is 

[ J d 3 xe lip '"- p ^ xa u(p', s'){y°(p> + q 3 /2) + E^} 

u(p,s)(V/2)} 
< jp xe *W a -P«)x a u 1[( J / } s '^ p 3 + q J/2)(V/2) 

+E(h j /2)a j }u(p,s)]. 

(7.16) 



Note that the term q ■ h will vanish by gauge choice, just 
as the q ■ A term vanished for the electromagnetic case. 
Then, using the same manipulations as previously on the 
spin products we reduce this to 

d 3 xe^- p « )xa u(p', s'HtV' + q 3 /2) + Ej j } 

u(p,s)(h j /2)] 

= [ J d 3 xe l ^~ p « )xa x\,{p- h+^\7xh- <t} X ], 

(7.17) 

where we have neglected terms of higher order, that is 
1/m 2 . Finally we combine (7.15) and (7.17) to obtain 
the total energy 

[J d s xe i ^- pa)x °'xl(m(l)+p-h+^7 x h-a) x ] 
= J d 3 x^(mcj) + p-h+^n- a)f. (7.18) 

(Recall that the gravitomagnetic field is O = V x ft.) 
This is the main result of this section and is consistent 
with the result of ref. [T5] for a scalar particle. 

Finally we note that since we have expanded the wave 
function in terms of a free Dirac particle on the mass 
shell (7.2) the free Dirac Lagrangian is zero and the extra 
geometric interaction term (h/2)L discussed in section 6 
vanishes. 



8. GRAVITOMAGNETIC PHYSICAL EFFECTS 

The result (7.18) is to be compared with the analogous 
electromagnetic result (7.10). We see, of course, that 
the Newtonian potential is the analog of the Coulomb 
potential eA° and the gravitomagnetic potential is the 
analog of the vector potential according to 



eA° o 0, (-e/m)A<+ h. 



(8.1) 



We also see that the coupling of the spin to the gravit- 
omagnetic field Cl is only half the analogous electromag- 
netic coupling. To make this most obvious we consider 
a gravitomagnetic field Q that is approximately constant 
over the system so that we may choose h = (f2 x r)/2. 
Then 



AE a = 



d 3 x^{mcf) + ~fi X f -p+ if2 • 
d 3 x^(m<j) + ^0 • r x p + ifi • ^)* 
d 3 xtf[m<t> + -ti- (L + S)}® (8.2) 
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Both orbital and spin angular momenta couple in the 
same way to the gravitomagnetic field, so there is no 
anomalous g factor for gravitomagnetism; that is g g = 1 
for both orbital and spin angular momenta. 

From the above correspondence it is clear that since 
a magnetic moment due to orbital angular momentum, 
(e/2m)L, precesses at the Larmor frequency (eB /2m) in 
a magnetic field B, the gravitomagnetic moment due to 
both orbital and spin angular momenta will precess in a 
gravitomagnetic field fl with frequency f2, but in the op- 
posite direction. Thus quantum precession should be the 
same as that observed in the classical gyroscope systems 
of the GPB satellite experiment [3]. It thus seems very 
likely that the precession rate is universal for any angu- 
lar momentum system, whether the angular momentum 
is classical or quantum mechanical, orbital or spin. 

For the surface of the earth the magnitude of the grav- 
itomagnetic field is quite small, as estimated in ref. [15] 
The field and the associated quantum energy are of order 



CI ps 10~ 16 rad/s, E n = hfl ps lO'^eV. 



(8.3) 



Experimental detection of such small quantum gravito- 
magnetic effects in an earth-based lab would obviously 
be difficult. Such an experiment might be performed 
with an atomic interferometer. The atomic beam could 
be split into two components with angular momenta 
differing by AL = h. Then, according to (8.2) the 
two components would have energies differing by about 
AE ps QAL pa Vlh and thus suffer phase shifts differing 
by about Aip ps AEt/h pa fit, where t is the time of 
flight. For a typical t = Is this implies a phase shift of 
order 10 _ 13 ra<i, which is orders of magnitude less than 
presently detectable [2"5] . 

In addition to the small size of gravitomagnetic effects 
one might see in the laboratory there is a serious fur- 
ther inherent difficulty in almost any such experiment; a 
rotation of the apparatus would in general have similar 
effects and swamp the gravitomagnetic effects, so such 
rotations would have to be controlled and compensated 
to very high accuracy as mentioned in the introduction 
and in ref. [15] . 

The results of the GPB experiment and the theoreti- 
cal results of this paper and ref. [15] are probably most 
important in establishing the validity and consistency of 
general relativity and the gravitomagnetic effects that it 
implies. Such gavitomagnetic effects are quite small in 
earth-based labs and satellite systems, as is clear from 
(8.3), but may play a large role in astrophysical phenom- 
ena such as the jets observed in active galactic nuclei, for 
which the gravitomagnetic fields are much stronger [S]. 



field and the frame-dragging or gravitomagnetic field 
take simple and intuitive forms. Detection of the small 
gravitomagnetic effects in earth-based or satellite exper- 
iments is quite difficult, but such effects are expected to 
be large and of great interest in astrophysical systems 
such as jets from active galactic nuclei and black holes. 



Appendix A. THE INVERSE DIFFERENTIAL 
OPERATOR 



We briefly study the type of differential operator that 
appears in (3.7) by solving the differential equation 

Af + df = (A + d)f = F, f = f(x), F = F(x), 

(A.l) 

where F(x) is a given function that may be expanded as a 
power series in the region of interest and A is a constant. 
The solution of the homogeneous equation is 

f h = Ce~ Ax (C = arbitrary constant). (A. 2) 

The general solution of (A.l) is fh plus any particular- 
solution f p ; for the particular solution we solve (A.l) 
symbolically as, 



f p = (A + d)- 1 F=±( 



1 / dF d 2 F 
A A 2 



(A.3) 



Operating on (A.3) with (A + d) obviously gives F. 

To further justify the above formal operations we may 
solve (A.l) in a different way. An integrating factor is 
easily seen to be e Ax , so 

d{e Ax f) = e Ax {A + d)f = e Ax F. (A A) 

Integration then gives the general solution 



/ 



-Ax 



-Ax' 



F(x')dx' + Ce 



-Ax 



(A.5) 



Since (A.l) is linear and F is assumed to be expandable 
in a power series we need only consider powers, F = x n . 
Then we easily evaluate (A.5) using integration by parts, 
to obtain 



J A \ A 



A 2 



n(n — l)x 



71-2 



1 



Ce 



-Ax 



(A.6) 



This agrees with the power series for f p given in (A.3). 



9. SUMMARY AND CONCLUSIONS 

We have developed the theory of a spin 1/2 Dirac 
particle in a Riemann space and its weak field limit in 
considerable detail. In the low velocity limit for the par- 
ticle the energies due to the Newtonian or gravitoelectric 



Appendix B. ENERGY MOMENTUM TENSOR 
FOR THE DIRAC FIELD 



We wish to obtain the energy momentum tensor for 
a Dirac field in flat space, which occurs in (6.8) and 



11 



(6.9) [H]. We begin with the Lagrangian (3.1) for the 
free Dirac field and work out the canonical energy mo- 
mentum tensor according to the Noether theorem; it is, 
up to a constant multiplier C , 

fIT fIT 

= C\a^%^ v -M>\ v ^\. (B.l) 

where we have omitted the term proportional to L since 
it is zero for a solution of the free Dirac equation. Using 
the fact that the Dirac and the Klein-Gordon equations 
are obeyed by tp we calculate the two divergences of this 
tensor to be 

T" u \fi = 0, T» v \v = C{b - a)[m 2 (^W) - i$"vf4>\v)] 

(B.2) 

If we choose b = a, as in the text, both divergences are 
zero and the tensor has symmetry in yb and vb. Moreover 
we may then consistently symmetrize T^ v and have 



To verify the result (B.5) we may calculate the divergence 
of to find, after some algebra, that it gives the correct 
Lorentz force, 



If 



(B.6) 



In the interaction Lagrangian (6.8) the energy momen- 
tum tensor is contracted with the symmetric h^ v so the 
symmetrization in (B.5) is not relevant. 



1 



(B.3) 



This has now been normalized so that in the low velocity 
limit 



■oo 



rmjjtp 



(B.4) 



Finally, to include the electromagnetic field we use the 
minimal substitution recipe id^ 



id^ - eA^ to get 



1 



(B.5) 



Acknowledgements 



This work was partially supported by NASA grant 8- 
39225 to Gravity Probe B and by NSC of Taiwan under 
Project No. NSC 97-2112-M-002-026-MY3. Pisin Chen 
thanks Taiwan's National Center for Theoretical Sciences 
for their support. Thanks go to Robert Wagoner, Francis 
Everitt, and Alex Silbergleit and other members of the 
Gravity Probe B theory group for useful discussions, and 
to Mark Kasevich of the Stanford physics department for 
interesting comments on atomic beam interferometry and 
equivalence principle experiments. Kung-Yi Su provided 
valuable help with the manuscript. Elisa Varani thanks 
Cavallo Pacific for encouragement and support. 



[1] R. J. Adler, "Gravity," chapter 2 of The New Physics 
for the Twenty-first Century, edited by Gordon Fraser, 
(Cambridge University Press, Cambridge UK, 2006). 

[2] C. Will, Theory and Experiment in Gravitational 
Physics, (Cambridge University Press, Cambridge UK, 
1981, revised edition 1993), see chapters 2 and 5. 

[3] H. C. Ohanian and R. Ruffini, Gravitation and Space- 
time, (W. W. Norton, New York, 1976), chapter 1. 

[4] C. W. F. Everitt et. al. Phys. Rev. Lett. 
106, 221101(2011); see also the GPB website: 
|http: / /einstein.stan ford.edu 

[5] C. Will, "Finally, Results from Grav ity Probe B" , at 
http : / /physics . aps. org / articles / v4/43 

[6] I. Ciufolini et. al. Class. Quantum Gravit. 14, 2701 
(1997); I. Ciufolini et. al. Science 279, 2100 (1998); I. 
Ciufolini and J. A. Wheeler, Gravitation and Inertia, 
(Princeton University Press, New Jersey, 1995), chap- 
ter 6; I. Ciufolini et. al., in General Relativity and John 
Archibald Wheeler, edited I. Ciufolini and R. A. Matzner, 



(Springer, Dordrecht, 2010), p. 371. 

[7] L. Iorio, H. I. M. Lichtenegger, and C. Corda, Astrophys. 
Space Sci., 331, 351 (2011), discusses the phenomenology 
of solar system tests of general relativity; I. Ciufolini and 
E. C. Pavlis, Nature, 43, 958, (2004). 

[8] K. S. Thorne, Near Zero: New Frontiers of Physics, 
edited by J. D. Fairbank, B. S. Deaver, C. W. F. Everitt 
and P. F. Michelson (W. H. Freeman, New York, 1988) 
p. 573; L. Stella and A. Possenti, Space Sci. Rev., 148 
(2009); R. D. Blandford and R. L. Znajek, Mo. Not. Roy. 
Astro. Soc, 179, 433 (1977); R. K. Williams, (1995, May 
15). Phys. Rev., 51(10), 5387-5427, (1995). 

[9] S. W. Hawking, Commun. Math. Phys. 18, 301 (1970). 
[10] N. D. Birrell and P. C. W. Davies, Quantum Fields in 
Curved Space, (Cambridge University Press, Cambridge 
UK, 1982). 

[11] R. J. Adler, P. Chen, and D. Santiago, Gen. Rel. and 

Grav. 33, (2001). 
[12] V. Nesvizhevdky et. al. Nature 415, 297-299, (2002), 



12 



[13] S. Dimopoulos, P. W. Graham, J. M. Hogan, and M. 
A. Kasevich, Phys. Rev. Lett. 98, 111102 (2007); S. Di- 
mopoulos, P. W. Graham, J. M. Hogan, M. A. Kasevich, 
and S. Rajendran, Phys. Rev. D 78, 12202 (2008). 

[14] S. Weinberg, Gravitation and Cosmology, (John Wiley, 
New York, 1972), see chapter 5 on effects of general rel- 
ativity on the motion of a classical particle. 

[15] R. J. Adler and P. Chen, Phys. Rev. D, 82 025004 (2010). 

[16] J. D. Bjorken and S. D. Drell, Relativistic Quantum Me- 
chanics, (McGraw Hill, New York, 1964), we use the con- 
ventions and notation for the Dirac equation in flat space 
given in chapter 3; J. D. Bjorken and S. D. Drell, Rela- 
tivistic Quantum Fields, (McGraw Hill, New York, 1965), 
chapter 13. 

[17] I. Lawrie, A Unified Grand Tour of Theoretical Physics, 
(Adam Hilger, Bristol, 1990), generally covariant Dirac 
theory using tetrads is discussed in sec. 7.5; see also ref. 
[14] sec. 12.5. 

[18] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravita- 
tion, (W. H. Freeman, San Francisco, 1970); R. J. Adler, 
M. Bazin, and M. M. Schiffer, Introduction to General 
Relativity, 2nd edition, (McGraw Hill, New York, 1975), 
see chapter 9. 

[19] R. J. Adler, Gen. Rel. and Grav. 31, 1837 (1999.) 
[20] R. J. Adler and A. S. Silbergleit Int. J. Th. Phys. 39, 1291 
(2000). Note that the definition of the gravitomagnetic 



field differs by a factor of 2 from our use in this work. 
[21] R. P. Feynman, Lectures on Gravitation, lecture notes by 
F. B. Morinigo and W. C. Wagner (Caltech, Pasadena, 
1971); See also ref. [H], sec. 6.9; P. C. Peters, Flat-space 
Gravitation and Feynman Quantization, lecture notes, 
(1965). 

[22] An overview of this vast field is given in D. Oriti, Ap- 
proaches to Quantum Gravity, (Cambridge University 

Press, Cambridge UK, 2009). 
[23] A. Trautman, "Einstein-Cartan Theory," in Encyclopedia 

of Mathematical Physics, edited by J. -P. Francoise, G. L. 

Naber, and S. T. Tsou, (Elsevier, Oxford UK, 2006), p. 

189; H. Kleinert, [gr-qc] arxiv: 1005. 1460 (2010) 
[24] M. E. Peskin and D. V. Schroeder, An Introduction to 

Quantum Field Theory, (Addison- Wesley, Reading MA, 

1995), chapter 3. 
[25] R. Shankar, Principles of Quantum Mechanics, 2 nd ed. 

(Plenum, New York, 1994), section 20.2. 
[26] J. Schwinger, Am. J. Phys. 42, 507 (1973); J. Schwinger, 

Gen. Rel. and Gravit. 7, 251 (1976). 
[27] R. J. Adler and S. D. Drell, Phys. Rev. Lett. 13, 349, 

(1964); R. J. Adler, Phys. Rev. 141, 1499-1508 (1966). 
[28] R. J. Adler, H. Mueller, and M. L. Perl, Int. J. of Mod. 

Phys. A, (2011); to be published. 
[29] See also chapter 5 in ref. [16] 



